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ABSTRACT 

We present stellar evolutionary sequences for asymptotic giant branch (AGB) 
stars in the Magellanic Cloud clusters NGC 1978, NGC 1846 and NGC419. The 
new stellar models for the three clusters match the observed effective tempera- 
tures on the giant branches, the oxygen-rich to carbon-rich transition luminosi- 
ties, and the AGB-tip luminosities. A major finding is that a large amount of 
convective overshoot (up to 3 pressure scale heights) is required at the base of the 
convective envelope during third dredge-up in order to get the correct oxygen-rich 
to carbon-rich transition luminosity. The stellar evolution sequences are used as 
input for detailed nucleosynthesis calculations. For NGC 1978 and NGC 1846 we 
compare our model results to the observationally derived abundances of carbon 
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and oxygen. We find that additional mixing processes (extra-mixing) are re- 
quired to explain the observed abundance patterns. For NGC 1846 we conclude 
that non-convective extra-mixing processes are required on both the RGB and 
the AGB, in agreement with previous studies. For NGC 1978 it is possible to ex- 
plain the C/0 and ^^C/^^C abundances of both the 0-rich and the C-rich AGB 
stars by assuming that the material in the intershell region contains high abun- 
dances of both C and O. This may occur during a thermal pulse when convective 
overshoot at the inner edge of the flash-driven convective pocket dredges C and 
O from the core to the intershell. For NGC 419 we provide our predicted model 
abundance values although there are currently no published observed abundance 
studies for the AGB stars in this cluster. 

Subject headings: Galaxies: star clusters: individual: NGC 1978, NGC 419, 
NGC 1846 - Magellanic Clouds - Stars: abundances - Stars: AGB and post- 
AGB - Stars: evolution - Nuclear reactions, nucleosynthesis, abundances 



Introduction 



Asymptotic Giant Branch (AGB) stars are evolved, low to intermediate mass stars (1 
< Mj/Mq < 8). They are characterized by a degenerate inert C-0 core surrounded by 
helium and hydrogen nuclear burning shells, burning alternatively, forming a double-shell 
configuration. An intershell region rich in helium and carbon exists between the He and H 
nuclear shells. A deep H-rich conve ctive envelope surrounds the double-shell configuration 
(e.g., see the review by iHerwi gl |2 00 5l ) . During the thermally-pulsing AGB (TP-AGB) phase, 
the He-shell becomes unstable, igniting every s^lO^ years or so, with the resulting thermal 
pulses (TPs) lasting for ?«10^ years. Following TPs mixing episodes can occur. The mixing 
brings the products of nuclear b urning of H and He (mostly ^He and ^^C) from the interior 
of the star to the stellar surface (Ilbenlll975l ). These are referred to as third dredge- up (TDU) 
events. Through the action of repeated TDUs, AGB stars evolve from having an oxygen- 
rich composition (0-rich stars) where C/0 < 1 to a carbon-rich composition (C-rich stars) 
where C/0 > 1. T he TDU can also r esult in AGB stars with enhanced s-process elements 
in their spectra (see IBusso et al.lll999h . Thus, thermal pulses lead to changes in the surface 
abundances of AGB stars, making them chemically very different compared to their less 
evolved counterparts. The end of the AGB phase is reached when the superwind mass loss 
(up to ~ lO~^M0yr~^) reduces the hydrogen-rich envelope to small value s (<1O~^M0), with 
the ejected matter diffusing into the inter-stellar medium jBlockeilboOll l 



Star clusters are ideal sites to test theories of stellar evolution. They contain stars 
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of similar age and metallicity. Star clusters in the Magellanic Clouds (MCs) prove very 
valuable in comparison to the star clusters in our Galaxy as they span a wide range of 



age, which enables us to study the evolution of stars of various ma sses (e.g. iGirardi et al. 



20091 : iMackev et all 120081 : Milone et al.ll2009l : iMucciareUi et al.ll2006l . etc). The MCs house 



a large population of rich, intermediate age clusters which are useful for studying the the 



(Girardi et al. 


1995; 


Wood 


1994) 



the luminosities of AGB stars can be determined, in contrast to AGB stars in the Galaxy 
which occur mostly in field popu lations. Further, the intermediate-ag e clusters in the MCs 



demonstrate the TDU brilliantly ( iBessell et al.lll983l : lFrogel et al.lll990[ ) and are useful probes 



to study the details of the TDU in 0-rich stars and C-rich stars. 

The objective of our work is to study the evolution and nucleosynthesis of AGB stars 
in the two Large Magellanic Cloud (LMC) clusters, NGC 1978 and NGC 1846, and in the 
Small Magellanic Cloud (SMC) cluster NGC 419. These three clusters are ideal test beds for 
evolution and nucleosynthesis studies owing to the availability of accurate estimates of AGB 
structural parameters such as pulsation masses, effective temperatures (Tefj), and luminosity 
(Kamath et al. 2010 for NGC 1978 and NGC 419, and Lebzelter & Wood 2007 for NGC 1846). 
Abundance studies and attempts to explain the observed C and O abun dances for NGC 1978 
and the C, O, and F abundances for NGC 1846 have been carried out by lLederer et al.l (|2009[ ) 
and iLebzelter et al.l (|2008[ ). respectively. It was found that the derived C/0 and ^^C/^^C 
ratios for the two clusters showed very different results considering that the two clusters have 
similar AGB masses and metallicities. For NGC 1978, no theoretical scheme was established 
that could satisfactorily reproduce the derived chemical abundance pattern. Further, it was 
found that the M-stars in NGC 1846 showed a v ery rapid increase in the observed [F/Fe] 
versus C/0 ratios compared to the predictions ( ILebzelter et al.ll2008l ). For NGC 419, no 
detailed observational information on abundances exist. 

In our study, we compute new evolution models with updated opacities for the AGB 
stars in these clusters. We constrain our models based on the accurate observational pa- 
rameters. This allows us to reproduce observables such as the giant branch temperatures, 
the oxygen to carbon transition luminosity (where C/0 ~ 1), and the AGB-tip luminosity. 
We then use these stellar evolutionary sequences in a post-processing code to study detailed 
nucleosynthesis and try to reproduce the observed abundances of the AGB stars in our target 
clusters. We estimate the effects of extra-mixi ng processes on the su r face abundanc e s, sirn - 
ilar to those atternpted i n previous studies by ILebzelter et al.l ( l2008l ). iLederer et al.l ( l2009l ). 
and lKarakas et al.l (120101 ). We also simulate the effects of an intershell enhanced in both ^^C 
and ^^O to explain the unusual abundances of the cluster AGB stars. 
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The outline of this paper is as follows. In Section [2] we introduce the three target clusters 
and supply existing abundance information for them. In Section [3] we present a general 
overview on the numerical method involved in the stellar evolution and nucleosynthesis 
models that we use in our work. Here we also discuss the uncertainties in stellar models and 
details on the updated stellar evolution code. In Section H] we present details of the individual 
models for each cluster and we discuss our results. Finally, in Section [5] we summarize our 
work and provide some concluding remarks. 



2. Target Clusters 
2.1. NGC1978 

The rich intermediate-age cluster NGC 1978 is an interesting candidate in the LMC 



This cluster houses both M-type and C-type stars (ILloyd Evand Il980l : iLederer et al.l 12009 



Kamath et al.ll2010[ ). From pulsation analysis , the cluster is found to have red variables early 
on the AGB with a mass of 1.55 ± 0.1 Mq (IKamath et al.ll2010l ). This study also showed 
that the highly evolved AGB stars have had a substantial ar nount of mass loss. NGC 1978 
is the only cluster in the LMC with a known mid-IR source (ITanabe et al.l 119981 ) which has 
a large infra-red excess, indicative of a very large mass-loss rate. 



The properties of the cluster are given in iKamath et al.l (120101 ). The metallicity esti- 
mates are mostl y in the range [Fe/H] = —0.37 dex to —0.42 dex with no evidence for any 
a-enhancement (iMucciarelli et al.ll2008l ) and the age is estimated to be r = 1.9 ± 0.1 Gyr. 
The age and metallicity estimates lead to an initial mass for the current AGB stars of 1.54 
to 1.62 Mq, consistent with the direct pulsation mass determinations for early- AGB stars of 
1.55 ± 0.1 M0. 



Lederer et al.l (120091 ) derived C/0 and ^^C/^^C ratios for nine AGB stars in this cluster. 



They found that, for the M-stars in their sample, the C/0 ratio ranges from 0.13 to 0.18 
with a typical uncertainty of ±0.05. The ^^C/^^C ratio values range between 9 and 16 with 
an uncertainty of up to ±4. They attributed these low C/0 and ^^C/^^C ratios for the 
M-stars in the cluster to the fact that the M-stars have not undergone any TDU. Based on 
a sample of four C-stars, for which reliable abundances could be established for only two, 
they found that the C/0 ratio was around 1.35 with an uncertainty of up to ±0.10 and the 
corresponding isotopic carbon ratio was about 150 to 175 with significantly larger errors (~ 
±25 to ±50). 
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2.2. NGC1846 



NGC 1846 is a rich interm ediate-age LMC cluster known to have a very interesting 
CMD. Using HST observations, iMackey &: Broby Nielsen! (j2007| ) identified the presence of 
two distinct main-sequence turn-offs that are clearly associated with the cluster. They are 
the result of the presence of two separate stellar populations with the same metallicities but 
different ages o f r„ = 1.5 Gyr for the u pper turn-off and ri = 1.8 Gyr for the lower turn-off. 
Recent work by lGoudfrooij et al.l ( l2009l ) using HST data confirmed the double main sequence 
turn-off feature as well as identified the presence of a RGB bump, a feature that is also found 
in NGC 1978. They found that age values of 1.7 Gyr and 2.0 Gyr for the upper and lower 
turn-off's respectively, AGB star masses between 1.70 and 1 . 77 M (t,, and best-fit abundances 
of [Fe/H] = -0.50 with [a/Fe] = 0.20. Ilebzelter fc Wood! J2007h derived pulsation masses 
for the AGB stars in NGC 1846 of I.8M0, which corresponds to a cluster age of 1.9 Gyr. 



Lloyd Evand (Il980l ) identified a large number of M-stars, C-stars, and S-stars in NGC 1846 . 



The cluster AGB stars do not show significant mass loss along the AGB (ILebzelter fc Wood _ 
2007). Further, no star s with a high mid-IR excess have been found in NGC 1846 (jTanabe et al 
1998f l. ILebzelter et~aD J2008h determined the C/0, ^^C/^^C and the [F/Fe] ratios for a small 
sample of AGB stars in this cluster. For the M stars, they determined C/0 ratios between 
0.2 and 0.65 with an uncertainty of up to 0.1 dex. Carbon isotopic ratios varying between 
12 and 60 were found for the sample of M-stars. For the C-stars in the sample, they derived 
a C/0 value of around 1.8 and an isotopic carbon ratio of about 60. For the M-stars in 
their sample, they also measured the change of fiuorine abundance along the AGB using the 
blended HF line and found a clear increase in the F abundance with luminosity with the 
[F/Fe] values ranging between —0.71 to 0.40. 



2.3. NGC 419 



NGC 419 is a populous intermediate-age cluster in the SMC. This clu s ter houses a large 



popu l ation of AGB stars, many of which are C-stars (iFrogel et al.l Il990 



Mucciare 



li et al. 



20081 : iKamath et al.ll2010l ). The properties of this cluster are given in lKamath et al.l ( 120101 ). 
The metallicity is estimated to be around [Fe/H] = — 0.7 dex and the ag e estimates are mostly 



in the range 1.2 to 1.6 Gyr. Using the variability of stars in this cluster, iKamath et al.l (120101 ) 
derived pulsation masses of 1.87 ± 0.1 Mq early on the AGB. The AGB stars in this cluster 
show significant mass loss along the AG B which agrees well with the existence of a mid-IR 
source detected by the ISO CAM survey (jTanabe et al.lll998r ). 
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The Numerical Method 



We calculate the stellar evolution and nucleosynthesis in two steps. First, we use the 
stellar evolution code to follow the evolution of the stellar structure and abundances impor- 
tant for stellar evolution (H, '^He, ^He, ^^C, ^^N, and ^^O) from the zero-aged main sequence 
(ZAMS) to the end of the TP-AGB phase (see Section [XT]) . Then we perform detailed nu- 
cleosynthesis calculations (see Section [3l2|) . The numerical method a nd the procedu r e use d 
to compute the mo d els have been p reviously described in detail by iKarakas et al.l ( 12002| ). 
Lugaro et al.l ( 120041 ). iKarakad (120101 ). and references therein. Here, we summarize the essen- 
tial details relevant to our study. 



3.1. Stellar Evolution models 



The stellar evolution is calculated using an updated version of the Mount Stromlo Stellar 

Evolu tion Code ( IWood fc Zarr 3l98l[lLattanziJl986l : IProst fc Lattanziolll996[|Karakas fc Lattanziol 
20071 ). The masses and compositions of the stars evolved are listed in Section |H Low- mass 
stellar models are affected by many uncertainties, the most important of which are listed 
below. With our improved models we aim to constrain these uncertainties. 



3.1.1. Mass-loss 



Dealing with the extent and temporal distribution of mass-loss in AGB stars is a ma- 
jor uncertainty in stellar modeling. Model calculations use simple parametrized formulae 



code, the 


Reimers 


1975) r 


Vassiliadis & Wood 


( 


1993) 



0.4, and the 



latter prescription was derived from a sample containing both 0-rich and C-rich stars so 
that it should be valid for our modeling. Furthermore, at higher mass-loss rates, direct 



measurements of the mass-loss rates for 



VIC AGB stars a grees reasonably well with the 



Vassiliadis fc WoodI (119931 ) mass-loss rate ( Wood et al.ll2007 ). In our study to reproduce the 
AGB-tip luminosity, we found that the IVassiliadis fc WoodI (119931 ) mass-loss prescription 
produced AGB-tip lumi nosities that were ab o ut 0.2 mag too faint. As a simple correction 
for this, we modify the IVassiliadis fc WoodI (119931 ) mass-loss prescription and invoke the 
super-wind phase at a later stage during the AGB evolution. This is done by effectively 
increasing the pulsation period at which the super-wind starts from 500 days (as found in 
Vassiliadis & Wood 1993) to ~ 710 — 790 days (specifically, the term 'P' in Equation 2 of 
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Vassiliadis fc WoodI ( 119931 ) is replaced by 'P — 210' and 'P — 290', respectively, depending 
on the mass and metallicity combination). The period values at which the superwind mass- 
l oss ra te begins (i.e., where the mass- loss rate given by Equation 1 of IVassiliadis fc Wood 
(119931 ) matches the value given by the modified Equation 2) for each structure model is listed 
in Table [H 

It is worth commenting on the need to increase the period for the onset of the superwind. 
The basic reason for this is that in these models with new C-rich opacities, the AGB stars 
become considerably cooler when C/0 exceeds unity. This means that the radius and hence 
computed pulsation period increases significantly causing the superwind to occur at lower 
luminosities than in 0-rich stars or in past models without C-rich opacities. Observatio nally, 
as noted above, the 0-rich and C-rich stars in the sample of IVassiliadis fc WoodI (119931 ) seem 
to fall on the same curve of rising mass loss with pulsation period. Also, it appears that 
the K — log P a nd Mhni — log P r elations are indistinguishable, at least while the stars are 



optically visible ( iFeast et al.lll989l ). There seems to be an inconsistency between the observed 
M — P relation and that predicted by the stellar models with the new C-rich opacities. Our 
adjustment of the onset period for the superwind is one way to bring the model mass loss 
rates and periods back to those that seem to apply observationally. Perhaps new pulsation 
models with C-rich opacities will resolve this problem. 



3.1.2. Convection 



One of the biggest uncertainties in stellar models is the treatment of convection. We 
use the mixing length (MLT) theory for convective regions and we set the mixing length 
parameter a = i/Hp in t he models by matchin g the observed RGB and E-AGB T^g, taking 
observed Tes values from iKamath et al.l (|2010( ) . We then keep the mixing length parameter 
a constant. 



Previous studi e s bvlStraniero et al.l ( 119971 ) , iKarakas et al.l (l2002h , IStancliffe et al.l (120041 ) , 
and Karakas et al. ( 201ol ) show that different stellar evolution codes predict differ ent TDU 
efficie ncies, with some codes predicting no TDU if overshoot is not included (e.g.. iMowlavi 
1999 ). Furthermore, low-mass (M < 2 Mq) AGB models show little or no TDU (IKarakas et al, 
2OO2I ) . In our models, to alter the extent of TDU, we include convective overshoot by ex- 
tending the position of the b ase of the convective envelope downward by A^ov pressure scale 
heights (jKarakas et al.ll2010l ). We include convective overshoot at the base of the enve- 
lope at all times during the AGB. By changing Aqv we alter the amount of carbon that is 
dredged- up to t he outer layers of the star such that the M/C transition takes place near the 
observed Mboi (ILebzelter fc WoodI l2007l : iKamath et al.ll2010l ). We note that we define the 
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M/C transition luminosity as the luminosity of the brightest M-star as there is an overlap 
in the luminosity of M and C-stars due to luminosity variations over a thermal pulse cycle. 
The values of a and A^^ov that we use for the models of the cluster AGB stars are listed in 
Table [H 

During thermal pulses, adding convective overshoot to the base of the intershell convec- 
tion will cause this convection to penetrate further into the C-0 co re, resulting i n intershell 



abundances that are different from models with no overshoot (e.g.. iHerwigl 120001 ) . The sub- 
class of PG1159 post-AGB stars are H-deficient and show He-intershell material at their 
surface (IWerner et al.ll2009[ ). In these stars, C abundances vary from 15-60% (by mass) and 
O from 2-20%. These abundances are in direct contrast to stand ard intershell compositions 
that giv e 25% C and < 2% O (e.g.. lBoothroyd fc Sackmannlll988l ). The theoretical models of 
Herwigl ( 120001 ) include the effect of diffusive convective overshoot into the C-0 core during a 
TP. This increases the C and O intershell abundances well above the value found in standard 
models but consistent with PC 1159 star abundances. Though we do not include this in our 
evolution models, in our nucleosynthesis study for AGB stars in NGC 1978 (refer Section l^?T]) . 
we synthetically estimate the effects of a non-standard intershell which is enha n ced in both 
^^C and ^^O . We employ methods similar to those described in iKarakas et al.l (120101 ). 



3.1.3. Opacities 



The outer layers of AGB stars become cool, allowing molecules to form. Therefore an 
accurate treatment of molecular opacities is needed. In our version of the stellar evolution 
code we have utilized new opacity tables. At low temperatures we ha ve use d the Ro s seland 
mean opacities computed using AESOPUS JMarigo fc Aringeil l2009h with lloddersl J2003h 
solar abundances as the reference solar composition. Using AESOPUS we have the option 
of incorporating scaled-solar, a-enhanced, and carbon-depleted opacity tables. We have 
generated tables based on the initial composition we require for ea c h of the clusters. We 
al so utilize the O PAL radiative opacity tables of llglesias fc Rogersl ( 119961 ) updated to use 



Lodderj ( l2003l ) solar abundance distribution of elements from C to Fe. To maintain 



consistency between the abundances of the low-temperature tables and the high-temperature 
OPAL tables we use OPAL tables appropriate for each of the compositions (scaled-solar, a- 
enhanced, and carbon-depleted). 
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3.2. The nucleosynthesis models 



The evolution code only includes the species (H, ^He, ^He, ^^C, ^^N, and ^^O) that are 
relevant to the major energy generating reactions and this forms the basis for the structure 
of the model star. In order to try and explain the elemental and isotopic abundance pat- 
terns observed in stars we need to include more nuclear species. We use a post-processing 
nucleosynthesis code for this purpose. This code needs as input from the stellar evolution 
code variables such as temperature, density, and convective boundaries as a function of time 
and mass fraction. The code then re-calculates the abundance chang function of 

mass and time using a nuclear network which contains 77 species (from hydrogen to sulphur, 
along with a small group of iron-p eak elements) and time dependant diffusive mixing for 
all convective zones (jCannonlll993l ). The overshoot regions at the bottom of the convective 
envelope in the evolution code are not treated as convective in the nucleosynthesis code, 
i.e., convection and mixing are assumed to stop at the Schwarzschild boundary. The code 
also requires input physics such as reaction rates and initial a bundances. Most o f the 589 
reaction rates are taken from the JINA REACLIB data bas e (ICvburt et al.ll2010l ). Details 
on the updated reaction rates that we use can be found in KarakasI ( 201ol ). For the three 
clusters, we assume the initial abundances based on the existing observational information 
on their abundances. Details on the input abundances u sed for each mo del are discussed in 
Section m We take our reference solar composition from iLoddersI (120031 ). 



An additional feature of the post-processing nucleosynthesis code which is not in the 
evolution code is that we can estimate the effect of extra-mixing. Observations of low- 
mass red gi ant stars (M < 2Mq near the RGB-tip) reveal ^^C/^^C ratios of ~ 10 and 
C/N ~ 1.0 (iGilroyl Il989l ). These ratios are lower than predicted by st andard stellar evo - 



lution models, which give ^^C/^^C ~ 20 and ^"^C/^^N ~ 1.5. (e.g. ICharbonnell Il994[ ). 
These trends are also seen in globular clusters whe re there is an anti-correlation of C 
abundance with luminosity (e.g., in M3; ISmithl 120021 ). Together these observations indi- 



cate the occurrence of non- convective mixing processes on the g iant branch. Mechanisms 



proposed to account for this extra-mixing fseelHerwig et al.ll2006f) in clude rotational mixing 



( lEggleton et al 



magnetic fields ( iNordhaus et al. 



20081: 



( Charbonnel et al. 1998), gravity waves (Denissenkov &: ToutI 2000 ). thermohaline m ixing 



Charbonriel fc Zahnf 



2008 



2007 



Stancliffe et al 



2009 



Busso et al.ll2007l : iPalmerini et al 



Stancliffell2010h . and 



20091 ). To simulate 



the effect of extra-mixing on the RGB, we take the envelope composition at the tip of the 
giant branch from the nucleosynthesis calculations and alter it such that the ^^C/^^C ratio 
equals the observed M-star ^^C/^^C ratio. This is done by decreasing the ^^C abundance, 
and increasing the abundances of ^^C and ^^N in the entire convective envelope according 
to changes expected for CN cycling, i.e., the total number of ^^C, ^^C, and ^^N nuclei are 
conserved. The modified models are subsequently evolved through core-helium burning and 
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on th e AGB. An extensive explanation on this procedure can be found in iKarakas et aL 



torn 



4. Models and Results 

In this section we present the stellar evolution and nucleosynthesis models for the AGB 
stars in our target clusters. The parameters used to construct the stellar models are listed 
in Table [TJ These include the initial mass and metallicity, the initial abundance pattern, the 
log Teff for the M-stars on the AGB at Mboi = —4.00 (log Tcff_4), the mixing-length parameter 
(a) used to fit the giant branch Teff, the mass on the early AGB (Me_agb), the amount of 
convective overshoot required to fit the M/C transition luminosity (Nov), the bolometric 
luminosity of the model star at the M/C transition (M^^*"), the pulsation period (P) where 
the superwind mass-loss rate begins, and the predicted bolometric luminosity at the tip of 
the AGB (Mb"sf-*'P). 

We use starting compositions with a variety of C/0 ratios in our evolution calculations. 
These C/0 ratios were chosen to broadly match the C/0 ratios observed for the cluster 
M-stars which lie on the early- AGB. For NGC419 we construct models using a scaled- 
solar abundance pattern as there are no observed C/0 ratios. However for NGC 1978 and 
NGC 1846 we experiment with three different initial compositions as there are observed C/0 
ratios: a scaled-solar composition; a carbon-depleted composition where [C/Fe] = —0.25 dex 
with scaled-solar values for the other elements (hereafter, the carbon depleted model); and 
an a-enhanced composition where [a/Fe] = +0.20 dex while for all the other elements we 
assume scaled-so l ar val ues (hereafter, the a-enhanced model). For all the above mixtures we 



use the iLodderd (120031 ) reference solar composition. 



In Table [2] we present some of the details of the stellar structure models. The first row 
for each cluster gives the initial mass and metallicity. Then we include for each model the 
number of TPs computed, the H-exhausted core mass (hereafter core mass) at which the 
TDU begins [M^^), the maximum TDU efficiency {X^^^, where A = AMd^edge/AMh, AMh 
being the amount by which the core mass has grown between the present and previous TPs), 
and the average value for the TDU efficiency (Aavg)- We also include the total amount of 
mass dredged into the envelope during the AGB lifetime (Mdredge), the maximum He-shell 
temperature (T^^^^^q) which is usually taken from the final TP, the core mass and the total 
mass at the final (the last computed) timestep (Mc(f) and Mtot(f), respectively), and the 
final interpulse period (rip(f)). 
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4.1. Evolution and Nucleosynthesis model results for NGC 1978 

Based on the estimated AGB pulsation mass of ~ 1.55 Mq for the cluster AGE stars (see 
Section 12. ip we construct evolutionary sequences starting from the ZAMS having an initial 
mass of 1.63 Mq with Z = 0.006 (similar to the observed global metallicity for this cluster, 
[Fe/H] ^ —0.4) and Y = 0.25. We construct stellar evolution models based on three initial 
compositions: scaled-solar, carbon-depleted, and a-enhanced. To match the composition of 
the lowest M-star C/0 ratio, we assume that the a-enhanced model also has a slight carbon 
depletion of [C/Fe] = —0.05 dex. Figure [1] shows the theoretical Hertzsprung-Russell (HR) 
diagram for the carbon-depleted model as an e xample. We also ove rplot the AGB variables 



and the non- variable red-giants in NGC 1978 ( iKamath et al.ll2010l ). We find that the evo- 
lutionary track is in good agreement with the observed star positions. The other models 
corresponding to the other two initial compositions have similar evolutionary tracks. The 
values for the mixing length parameter [a) are the same for all three compositions at 1.90. 
However the amount of overshoot varies slightly (A^^ov ~ 2.45 — 3.00), with the a-enhanced 
model requiring the most overshoot and the carbon-depleted model requiring slightly less 
overshoot than the scaled-solar model (refer to Table [1]). The amount of overshoot used in 
each model series results in predicted M/C transition Mboi values th at give a good matc h 



to the observed M /C transition bolomet ric luminosity of Mbd = —4.5 (IKamath et al.ll2010l ). 



Observations from IKamath et al.l (120101 ) also indicate that the most luminous cluster AGB 
star has an Mboi = —5.06. To match this luminosity, the super-wind phase for the three 
models start at a pulsation period of 790 days. 

One interesting feature of NGC 1978's CMD is the presence of the RGB bump. This is 
because the RGB bump luminosity can be used as an extra observational constraint on the 
stellar structure models. The RGB bump is caused by the H-shell erasing the abundance 
discontinuity left by the retreating convective envelope during the first dredge- up (FDU). In 
Fig. [2] we show the location (in mass) of the inner edge of the convective envelope during 
the FDU. The black solid-line shows the scaled-solar model, and this model has a slightly 
higher RGB bump bolometric luminosity of Mboi = —0.53 when compared to the observed 
Mhoi ~ —0.31, calculat ed from the position of the RGB bump at V555 = 19.10 ± 0.10 



( iMucciarelli et al.ll2007f ). To match the observed RGB bump luminosity we can apply con- 
vective overshoot in a similar way as done for the AGB. The red dashed-line in Fig. |2] shows 
the result of one such model with a = 0.30. The figure shows that the depth of the 

FDU is not significantly altered, although the FDU begins at an earlier time and this lowers 
the RGB bump bolometric luminosity to better match the observed values. The structural 
details (e.g., surface luminosity, core mass, effective temperature, etc) were essentially the 
same at the tip of the RGB and on the E-AGB in the model with slight overshoot compared 
to the model without. For this reason we will ignore overshoot on the RGB for the rest of 
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Fig. 1. — The HR diagram for the carbon-depleted model for NGC1978. The cyan/grey 
lines mark the M/C transition and indicates that the star is C-rich. The black open circles 
represent the observed positions of the M-stars in NGC 1978, the black filled circles represent 
the observed positions of the C-stars i n NGC 1978, and th e smaller black dots represent the 
non- variable red-giants in the cluster (IKamath et al.ll2010l ). Note that the low luminosity of 
the C-star indicates that the star is in a post-flash luminosity dip. 
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this study. 

From Table |2] we find that all the three models in NGC 1978 experience ~ 15 TPs. 
To reproduce the observed M/C transition luminosity we require very extended TDU with 
A/'ov ~ 2.45 - 3.00. This resuhs in A^ax of 0.72 - 0.82 (see Fig. ED and an average A of 0.58 
— 0.65 for all TPs. When convergence difficulties terminated the evolution calculations (at 
an envelope mass of the order of ~ 0.05 to O.2OM0 for the different models computed), the 
mass-loss rates were so high that we would not expect any further TPs or TDU events to 
occur on the AGB. 

We now discuss the results of the post-processing nucleosjTithesis calculations. In Ta- 
ble [3] we include C, O, and F surface abundance information for all the stellar models. We 
provide the initial, post-FDU, and final computed C/0 ratios (C/Oi, C/Opost-FDU; and 
C/Of, respectively), where all abundance ratios are number fractions. Furthermore we list 
the initial, post-FDU, and final computed ^^C/^^C ratios (^^C/^^Ci, ^^C/^^Cpost-FDU, and 
i2q^i3q^^ respectively). In the last column we provide the [F/Fe] ratio at the tip of the AGB 
where we use the standard notation [X/Y] = log;^o(-^/Y)star ~ logiol-^/Y)©- 

From Table |3] we see that the scaled-solar model gives a post-FDU C/0 = 0.33, signif- 
icantly higher than the highest M-star C/0 ratio of 0.18. On the other hand the carbon- 
depleted models all yield a C/0 ratio of ~ 0.18, within the range of the M-star observations 
(C/0 = 0.13 to 0.18). However, the predicted post-FDU ^^C/^^C ratio does not match the 
observed range of ^^C/^^C = 9 to 16 for all carbon depleted models. In order to obtain a 
match, we assume extra- mixing has taken place on the RGB (as explained in Section I3l2|) . 
We alter the envelope composition at the tip of the RGB such that ^^C/^^C = 13, which 
agrees well with the observed average M-star ^^C/^^C value of ~ 13. We designate AGB 
evolutionary sequences in which the abundance change has been made as extra-mixing se- 
quences. 

In Fig.Hlwe show the evolution of the C/0 versus ^^C/^^C ratios for the carbon-depleted 
models on a logarithmic scale. I ncluded on the plot are the observed values for the M and 



C-type stars in NGC 1978 from iLederer et al.l (|2009[ ). The black solid-line represents the 



carbon-depleted model without extra-mixing on the RGB and the red dashed-line represents 
the carbon-depleted model with extra-mixing on the RGB. Furthermore, we also perform 
one nucleosynthesis calculation with a lower init ial ^^C/^^C = 50 c ompared to the solar value 



of 89 (from LMC HII regions, see discussion in lWang et al.ll2009l ). From Fig. IHwe see that 
both the models with no extra-mixing, i.e., either with a solar ^^C/^^C (black solid-line) 
or with an initial ^^C/^^C = 50 (blue dotted- line) , have a post-FDU abundance of ^^C/^^C 
which is higher than the observed M-star value. However, the model with extra-mixing (red 
dashed-line) has a post-FDU abundance of ^^C/^^C = 12.8 which by construct agrees well 
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Fig. 2. — Mass interior to the inner edge of the convective envelope, showing the depth of 
FDU for the scaled-solar NGC 1978 model. The black solid line represents the standard model 
without overshoot, and the red dashed-line indicates the case with mild RGB overshoot. 




Fig. 3. — Core masses for the scaled-solar NGC 1978 model. The black solid line shows the 
variation of the H-exhausted core mass with time, and the red dashed-line indicates the time 
variation of the He-exhausted core mass. 
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with the observed average M-star C/ C value of ~ 13. 

We also experiment with two a-enhanced models, one without extra-mixing on the RGB 
and one with extra- mixing on the RGB (see Table [3]), which result in post-FDU C/0 and 
i2q^i3q j-a^iQg similar to the c arbon-depleted case. Ho wever, NGC 1978 shows no evidence 



for an initial a-enhancement (IMucciarelli et al.l l2008l ) . Thus, we favour a model with a 



carbon-depleted initial composition coupled with extra-mixing on the RGB since it accurately 
reproduces the observed average C and O composition of the cluster M-stars. Howev e r, this 



scenario does not adequately explain the cluster C-stars. Previously, iLederer et al.l (120091 ) 
attempted to reproduce the observed carbon and oxygen abundances for the cluster AGB 
stars using models with an E-AGB mass = 1.55M0 and Z = 0.006 along with moderate 
extra-mixing on the RGB coupled with an [0/Fe] = 0.2 dex. These models were able to 
match the observed average C/0 and ^^C/^^C for the M-stars but they were unable to 
reproduce C/0 and ^^C/^'^C for the C-stars. 

The C-stars in NGC 1978 have large observed carbon isotopic ratios indicating a strong 
contribution of ^^C from the He-intershell mixed into the envelope by TDU. We note that the 
models show a slope very close to 1 implying that essentially only ^^C increases with time: 
any substantial change in slope would require an additional change to either ^'^C or O. One 
possible way to do this would be to have an intershell abundance with more ^®0 than that 
predicted by standard models as discussed in Section 13.1.21 In Figure [5] we show a synthetic 
AGB model prediction (red dashed-line) that results from increasing the C content to 40% 
and the O content of the intershell to 15% (by mass). The standard values for this model 
are 15% for C and 0.35% for O (by mass). The artificially ei ihanced numbers are wit hin the 



range of C and O abundances observed for PC 1159 stars (jWerner fc Herwigi 120061 ). With 
the enhanced intershell C and O values, the total elemental C and O surface abundances 
increase as a result of such mixing by factors of ~177 and ~10, respectively, in comparison 
to a model with a standard intershell which results in an increase by factors of ~ 86 and 
~1.05 in the surface composition. From Figure [5] we see that in NGC 1978 such a model 
fits the data for both the M-stars and the C-stars very well. Note that the final C/0 ratio 
(C/0 ~3.00) is much lower than the final C/0 ratio of the model with a standard intershell 
composition (C/0 ~ 7.60). This implies that some AGB stars may have intershells with 
higher abundances of ^^0 and ^^C than those predicted by standard models. 



4.2. Evolution and Nucleosynthesis model results for NGC 1846 



For the AGB stars in NGC 1846 we construct evolutionary sequences starting from the 
ZAMS having an initial mass of I.86M0 with Z = 0.006 and Y = 0.25. This results in a pre- 
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Fig. 4. — Log C/0 versus log ^^C/^^C ratios for NGC 1978's c arbon-depleted mod eL Also 
shown are the observational data for AGB stars in NGC 1978 ( iLederer et al.ll2009l ). where 
black filled triangles represent the M stars and the green filled squares represent the C 
stars. The black solid-line refers to the model without any extra-mixing, the red dashed-line 
represents the model with extra-mixing on the RGB, and the blue dotted-line has an initial 
i2q^i3q _ rjij-^g crosses in this plot and subsequent plots denote the envelope abundances 
after each TDU episode. 
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Fig. 5. — Same as Fig. IHbut showing results for NGC 1978's carbon-depleted model with 
extra- mixing (black solid- line), and the synthetic AGB model with an intershell enhanced 
in carbon and oxygen (red dashed-line). In the synthetic model we have assumed intershell 
mass fractions of 15% for O and 40% for C. 
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dieted E-AGB mass of 1.8 0M(7^ whieh is similar to the estimated pulsation mass of ~1.8OM0 
for the cluster AGB stars (ILebzelter &: Woodll2007l ). We compute stellar evolution sequences 
with three initial compositions: scaled-solar, carbon-depleted, and a-enhanced. Figure M 
shows the theoretical HR diagram for the g-en hanced model as well as the AGB variables 
and the non- variable red-giants in the cluster (ILebzelter fc WoodI 120071 ). We find that the 
observations generally match the theoretical evolutionary track well although the most lumi- 
nous M-stars are bluer than predicted by the evolutionary tracks. The HR diagrams for the 
models corresponding to the other two abundances mixes are similar. In Tabled] we list the 
input parameters required for the three sequences. We find that the value for a is the same 
for all three compositions at 1.74. Overshoot parameters of N^y ~ 1.05 are needed to re pro- 
duce the observed M/C transition luminosity of Mboi = —4.78 (ILebzelter fc Woodll2007l ) for 
the scaled-solar and the carbon-depleted model while the a-enhanced model requires slightly 
more overshoot (N^y ^ 1 .41). The AGB-tip bolometric luminosity for NGC 1846 is ~ —5.18 



(ILebzelter fc WoodI 120071 ). To match this luminosity the superwind phase had to start at a 
pulsation period of ~ 710 days. 

The models for NGC 1846 experience around 16 TPs (see Table [2]). The NGC 1846 
models have slightly less efficient TDU (Aavg ~ 0.46) when the observed M/C transition 
luminosity is reproduced, when compared to the NGC 1978 models (Aavg ~ 0.62) consistent 
with the smaller amount of overshoot (A^^ov ~ 1.05 — 1.41 compared to ~ 2.54 — 3.00 in 
NGC 1978). For the NGC 1846 model sequences, the final envelope masses when convergence 
difficulties terminate the evolution calculations range between ^ 0.075 to O.2OM0. We would 
not expect any further TPs or TDU episodes. 



Observations indicate that NGC 1846 shows a mild a-enhancement (IMucciarelli et al. 



20081 ) . Thus we choose the a-enhanced model as our favored model. Post-processing nucle- 
osynthesis calculations for the a-enhanced model start with an initial surface value of C/0 
= 0.31 and = 89, and yield post-FDU values of C/0 = 0.19 and ^^C/^^C = 22. 

We then alter the envelope ^^C/^^C composition such that ^^C/^'^C = 14 at the tip of the 
RGB, in a similar manner as for NGC 1978 (Table [3]). The M-stars in the cluster show a 
spread in C/0 ratios ranging from 0.20 to 0.65 and ^^C/^^C ratios between 12 and 60. In 
Figure [7] we show results from the a-enhanced models with extra-mixing on the RGB (red 
dotted-line) and without (black solid-line). From Figure [7] it is clear that an a-enhancement 
combined with extra-mixing on the RGB gives a good match to both the C/0 and ^^C/^^C 
ratios of the cluster M-stars. However, the subsequent AGB evolution from the M-stars does 
not explain the observed abundances for the cluster's C-stars. 



Stellar evolution calculations by ILebzelter et al.l (120081 ) show that by assuming extra- 
mixing on the RGB along with an initial [0/Fe] = 0.20 dex, the predicted carbon and 
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Fig. 6. — The theoretical HR diagram for the a-enhanced model for NGC 1846. The 
cyan/grey lines indicate that the star is C-rich. The black open circles represent the ob- 
served positions of the M-stars in NGC 1846, the black filled circles represent the observed 
positions of the C-stars in the NGC 1846, and the smal ler black dots represent the non- 
variable red-giants in the cluster (ILebzelter fc WoodI 120071 ). Note that the low luminosity of 
the C-star indicates that the star is in a post-flash luminosity dip. 
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Fig. 7. — Log C/0 versus log ^^C/^'^C ratios for NGC 18 46 's g-enhanced mode l. Also shown 



are the observational data for AGB stars in NGC 1846 (ILebzelter et al.l 120081 ) . where black 
filled triangles represent the M-stars and the green filled squares represent the C-stars. The 
black solid-line refers to the a-enhanced model without any extra-mixing, the red dashed-line 
represents the a-enhanced model with extra-mixing on the RGB. 
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Fig. 8. — The same as Fig. [7] for NGC 1846's additional model sequences. The black solid- 
line represents the scaled-solar model without extra-mixing on the RGB. The red dashed-line 
represents the carbon-depleted model without extra-mixing on the RGB and the blue dotted- 
line represents the carbon-depleted model with extra-mixing on the RGB. 
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oxygen abundances match the observed M-star comp ositions. This r esult agrees well with 



our conclusions. Similar results were also presented by lKarakas et al.l ( 120101 ) for the cluster's 



M-sta rs. In order to match the carbon and oxygen abundances in the C-stars, iLebzelter et al. 



( 120081 ) estimated the effects of extra-mixing on the RGB as well as a moderate extra-mixing 
in the late part of the AGB when the stars become C-rich. The theoretical models match 
both the observed M and C-star compositions as a result of the extra ^^C brought to the 
surface by extra- mixing on the AGB. However it is unclear as to why extra- mixing starts 
on the AGB only when C/0 > 1. Nevertheless we also conclude that extra-mixing may be 
required on the AGB to explain the compositions of the C-stars. It is possible that this extra- 
mixing on the AGB is weak hot-bottom burning where the base of the convective envelope 
extends into the upper regions of the H-shell and some CN-cycling occurs. Since hot-bottom 
burning increases with stellar mass, and the AGB stars in NGC 1846 are more massive than 
the AGB stars in NGC 1978, this could explain why this effect is seen in NGC 1846 and not 
NGC 1978. 

In Figure E] we show the results for the log C/0 versus log ^^C/^^C ratios for several 
additional model sequences. The scaled-solar model (black solid-line) does not match the 
C/0 and the ^^C/^^C ratios of the M-stars with the lowest C/0 ratios. These stars show 
no C-enrichment and probably have a composition similar to (or that of) stars at the tip of 
the RGB. We find that the scaled-solar model matches the C/0 and the ^^C/^'^C ratios of 
the other M-stars that show an enrichment of carbon caused by TDU, although their C/0 
< 1. However, we find that this model does not reproduce the ^^C/^'^C of the C-stars. For 
the carbon-depleted model without extra-mixing (red dashed-line) we find that though the 
predicted C/0 ratios lie in the range of the observed C/0 ratios for the M and the C-stars, 
we cannot match the ^^C/^^C ratio of the M-stars, where the post-FDU ^^C/^^C ratio for 
this model is 22. Furthermore, the ^^C/^^C ratio predicted for the C-stars is very far from 
the observed value. The carbon-depleted model with extra- mixing fits the observed C/0 
and the ^^C/^^C ratios of the M-stars in the cluster (with a predicted post-FDU ^^C/^^C = 
12, refer Table |3]). However, as for the previous models, this model also fails to reproduce 
the observational ^^C/^^C ratio of the cluster's C-stars. 



4-2.1. Fluorine Abundances in NGC 1846 



An additional observat ional constraint i n NG C 1846 is the [F/Fe] ratio (ILebzelter et al. 



20081 ). Observations from ILebzelter et al.l (|2008[ ) show a steep increase in the estimated 



fluorine abundance with the C/0 ratio. The black solid-line in Figure [9] shows the predicted 
[F/Fe] versus C/0 ratios using our best fit model (the a-enhanced model with extra-mixing) 
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for the AGB stars in NGC 1846. We also compute another model using the lowest observed 
[F/Fe] abundances (ILebzelter et al.l l2008l ) as a starting point, where the initial [F/Fe] = 
—0.71 (denoted using red dotted-hnes in Figure [H]). We find that the increase in the predicted 
abundance of fiuorine with the C/0 ratio from both the theoretical models is shallower than 
the increase in the obser ved fiuorine abundance for a given C/0 ratio. Theoretical models 
by ILebzelter et al.l (120081 ) also show a similar trend in the [F/Fe] versus C/0 ratios. 



Flourine production in AGB stars is quite complicated: F is synthesized during thermal 
pulses via a combination of neutron, proton, and alpha capture reactions. The most likely 
path for the product i on of fiuorine is via ^ ^N(a,7) ^^F(/3+)^^0(p,a)^^N(Q;, 7)^^F reactions 
( iForestini et al.l Il992l : iMowlavi et al.lll996[ ). The required ^^N can be synthesized by the 
^®0(p,a)^^N reaction owing to the presence of ^^O and protons in the He intershell. ^'^O 
is produced by the a-capture on ^"^N where ^^N is found in the ashes of CNO cycling from 
the preceding H-burning stage. Further, the ^^N(n,p)^^C reaction has a high cross section 
and can produce ^'^C and free protons. ^^F production is enhanced by the inclusion of a ^^C 
pocket via the release of free neutrons which come from the ^^C(a,n)^^0 reaction. At the end 
of each TDU episode, the convective envelope penetrates into the stable radiative intershell 
zone and a ^^C-rich region can form in the top layers of the He- intershell as a result of the 
partial mixing of protons. The neutrons from the ^^C(a,n)^^0 reaction not only produce the 
s-process elements but also get captured by species such as ^^N and ^^Al; both of which are 
strong neutron absorbers and are produced by the H shell. Neut ron captures on ^^N then 
help increase the ^^F production during a thermal pulse (see, e.g.. lLugaro et al.l 12004 ). 



In our best-fit model for the AGB stars in NGC 1846, we artificially include a partial 
mixing zone at deepest extent of e ach TDU episode in the nucleosynthesis code using the 
same procedure as that explained in Lugaro et al.l (2004) and Karakas ( 2010 ). which is based 
on the method used by lGoriely fc Mowlavil ( l2000l ). We insert the partial mixing zone at the 
deepest extent of the TDU because this is when a sharp discontinuity is produced between 
the convective envelope and the radiative intershell, which is a favourable condition for the 
occurrence of mixing. The partial mixing zone can be defined in terms of the proton profile. 
It is defined as the region where the abundance of protons drops exponentially from the 
envelope value to a fixed lower value (Xp). We experiment with two values of Xp = lxlO~^ 
and Xp = lxlO~^. For the size of the partial mixing zone we consider three values of 
MpMz = 4xlO-='Mo (i.e. 
and MpMz = exlO^^Mg 
TDU), and MpMz = 1.2xlO~^M0 (i.e., most of the mass of the intershell during the last 
few TP's). Figure [TU] shows MpMz, Xp, and the proton profile in two cases. Note that MpMz 
stays constant with evolution along the AGB. 



^1/6 of the mass of the intershell at the deepest extent of TDU) 
i^i.e., ~l/4 of the mass of the intershell at the deepest extent of 
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Figure [TT] shows the predicted [F/Fe] versus C/0 ratios for our best fit model for 
NGC 1846 after the inclusion of partial mixing zones starting with an initial [F/Fe] ~ —0.71. 
The red dashed-line represents the model where the partial mixing zone has a depth of MpMz 
= 4xlO~^M0 and a proton abundance limit, Xp = 1x10"'*. This model does not reproduce 
the observations. We experiment by extending the depth of the partial mixing zone to MpMz 
= 6xlO~^M0 and Xp = lxlO~^ denoted by the blue dotted-line in Figure [TH The model 
shows that a little more F is produced at a given C/0 ratio than the previous case, however, 
the fluorine abundance does not match the observations. The black solid-line in Figure [TT] 
represents the best-fit model of NGC 1846 with a Mpmz = 1.2x1O"2M0 and Xp = 1x10"^ 
This model reproduces the observed F abundances fairly well. This implies that a MpMz 
= 1.2xlO~^M0 or higher would be required to reproduce the observations: note that the 
total model intershell mass is reduced to ~1.5xlO~^M0 at the tip of the AGB. However, 
this situation is rather speculative as s-process models require smaller ^^C pockets in order 
to match observations (IGallino et al.lll998l : iBonacic Marinovic et al.ll2007l ). 

We conclude that our models can reproduce the observations but only with partial 
mixing zones that are likely to be larger than those required for s-process studies, making 
this a speculative result. We also note that the observed abundances were estimated based 
on the single blended HF line for the M-stars which also adds to the uncertainity in the 
observed F abundance. 



4.3. Evolution and Nucleosynthesis model results for NGC 419 



The last cluster in our study is the SMC cluster NGC 419. In Figure [T^ we show the the- 
oretical evolutionary track of the scaled-solar model and the posit ions of the AGB variables 
and the non- variable red-giants in NGC 419 (IKamath et al.l 120101 ). The AGB evolutionary 
track is clearly in good agreement with the observed star positions. Here we only consider 
a scaled solar initial composition owing to the lack of information on the abundances for 
stars in this cluster. We use a ZAMS model with an initial mass of I.9IM0, Z = 0.004, and 
Y = 0.25. This results in an E-AGB mass of I.85M0 which agrees well w ith the estimated 
pulsation mass for the AGB variables in this cluster (IKamath et al.ll2010l ). We find that an 
a of 1.74 is required to reproduce the giant branch temperatures: this is similar to that for 
NGC 1846 which is of a similar initial mass (Table [T]). An overshoot of N^v = 2.10 is re- 
quired on the AGB to reproduce the observed M/C transition luminosity (Mboi = —4.5). To 
match the observed AGB-tip luminosity, the super-wind mass-loss rate needs to begin at a 
pulsation period of 790 days. Owing to a lower metallicity and a higher mass, the AGB stars 
experience more TPs (19) compared to the two other cluster AGB stars (~ 15 each). The 
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Fig. 9. — [F/Fe] versus C/0 for NGC1846. The black solid line corresponds to the a- 
enhanced model with a initial scaled-solar [F/Fe] = 0.0. The red dotted-line corresponds to 
the a-enhanced model for which the initial [F/Fe] = —0.71. The b lack triangles denote the 
observed [F/Fe] values varying with C/0 for the cluster's M-stars ( iLebzelter et al.ll2008l ). 
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Fig. 10. — Proton profiles for tlie NGC 1846 a-enlianced model witli tlie inclusion of partial 
mixing zones after the 4**^ thermal pulse. The blue dashed-line shows the proton profile at 
the base of the convective envelope (shaded region) for MpMz = Qx10~^Mq, where Xp = 
1 X 10^^. The black dotted-line shows schematically the lower edge of the partial mixing zone 
in a case where MpMz = 4xlO~^M0, where Xp = 1x10^^. Also shown are the composition 
profiles for ^^C, ^^C, and ^^F. Note: The increse in ^'^C at ^ 0.586 is due to the start of the 
formation of the ^^C pocket. 
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Fig. 11. — The effect of partial mixing zones (pmz) and proton-profiles (Xp) on the [F/Fe] 
versus C/0 for NGC 1846 a-enhanced model. The black solid-line represents the model with 
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greater number of TPs leads to a large amount of He-shell burning material being mixed into 
the envelope and this is demonstrated by the high final C/0 and ^"^C/^^C ratios (Table |3]). 
The envelope mass when convergence failed is ~ O.29M0 which implies that the models star 
is unlikely to experience any further TPs or TDU episodes since the mass-loss rate at this 
stage is 9.4OxlO~®M0/year, and the interpulse period is ~ 1x10^ years. 



4.4. The effect of C/O ratio on T^s of AGB stars 



Figure dSl and Figure [H] show the variation of Mboi with log T^g measured at the 
interpulse luminosity maximum for each of the TP's for the two LMC clusters NGC 1978 
and NGC 1846 (black solid-lines). The size of the ellipses on the solid lines denotes the C/O 
ratio after each TDU episode, as indicated in the grid on the bottom-right corner of the 
figures. Also shown are the AGB variables that have observational C/O estimates (Lederer 
et al. 2009 for NGC 1978, and Lebzelter et al. 2008 for NGC 1846), with the M-stars depicted 
as black filled triangles and C-stars as green filled squares. The size of the ellipse around each 
of the filled symbols indicates the observational C/O ratio of the star. The AGB variables 
without any abundance information are denoted as open symbols (M-stars as black open 
triangles and C-stars as green open squares). For rn ost of the observati onal data, the Mhni 
and t he log T^g is taken from iKamath et al. J2010h for NGC 1978 and Ihebzelter fc Wood I 
(120071 ) for NGC 1846 and for those stars not in th ese studies, the M^and the log T^s is 
taken from lhederer et all ^00^) for NGC 1978 and Ihebzelter et al.l J2008h for NGC 1846. 



Figures [T3l and [T^ shows that for the M stars with luminosities —4 < Mboi < —3, T^s of 
the models matches the observed values of T^q well. This is because the mixing length was 
adjusted to provide a match between model and observed T^s values at these luminosities. 
At more luminous Mboi, the model Teg is slightly cooler than the observations. This can be 
attributed t o the effect of the consta nt mixing length that we use in our models. Previous 
studies (e.g. lLebzelter &: Woodll2007l ) have shown that at high luminosities the mixing length 



needs to be increased slowly with luminosity in order to reproduce the correct slope for the 
theoretical giant branch. 

Once the M/C transition luminosity is reached, it can be seen that the AGB slope of 
the models suddenly decreases and the stars become increasingly cooler with increasing C/O 
ratio. The few C-stars with observed C/O ratios generally show the same movement to 
lower Teff values as seen in the models, although the correspondence between Tes and C/O 
ratio is not clear. The larger sample of C stars that do not have C/O measurements also 
lie cooler than the sequence of 0-rich AGB stars or its extrapolation to higher luminosities. 
However, these stars do not necessarily lie directly on the sequences shown in Figures and 
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Fig. 12. — Evolutionary track for NGC419's scaled-solar model. The cyan/grey lines indi- 
cates that the star is C-rich. The black open circles represent the observed positions of the 
M-stars in NGC419, the black filled circles represent the observed positions of the C-stars 
i n NGC419, and the small black dots represent the non-variable red-giants in the cluster 
(IKamath et al.l 120101 ). Note that the low luminosity of the C-star indicates that the star is 
currently in a post-flash luminosity dip. 
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for models at the interpulse luminosity maximum because of the variation of luminosity 
over the thermal pulse cycle. This is best seen in Figures [H |6] and [121 In summary, the 
increase in C/0 ratio above unity causes T^s to decrease but it is hard to determine from 
the current observational data whether the magnitude of this effect is reproduced accurately 
by the models. 



Discussion and Conclusions 



In this paper we have presented new stellar models for AGB stars in NGC 1978, NGC 1846, 
and NGC 419. These are three clusters whose well-defined observational characteristics pro- 
vide strong constraints on evolution models. The stellar evolution models are constrained to 
reflect the observed AGB pulsation mass, cluster metallicity, giant branch effective temper- 
ature, M/C transition luminosity and the AGB-tip luminosity. 

Stellar models are affected by major uncertainties including mass loss, convection (e.g., 
mixing length), depth and efficiency of TDU, and non-convective mixing. In our study we try 
to constrain these u r icerta inties. In order to produce the correct AGB-tip luminosities, the 
Vassiliadis fc WoodI ( 119931 ) mass-loss prescription had to be modified so that the super-wind 
mass-loss rat e commences at a per i od of ~ 710 — 790 days as opposed to the 500 days period 
employed by IVassiliadis fc WoodI (Il993[ ). The mixing length parameter (a) required to fit 
the observed giant branch temperature for each of the three clusters shows values similar to 
the a's found when calculating a standard solar model (e.g.. iPiersanti et al.l 120071 although 
our values are slightly lower than their a ~ 2.1). 

Overshoot at the base of the convective envelope of ~ 1 — 3 pressure scale heights was 
required to obtain the observed luminosity of the M/C transition. This results in higher final 
C/0 ratios of about ~ 7.5 - 8.0 for the LMC cluster NGC 1978, ~ 3.7 - 4.8 for the LMC 
cluster NGC 1846, and ~ 11.9 for the SMC cluster NGC419. The C/0 ratios of planetary 
nebulae (PNe) pu t uppe r limits on the C/0 ratios for the LMC and SMC. In the LMC, 
Stanghellini et al.l ( 120051 ) found that a C/0 of 2 is typical for non-bipola r PN and in the 
SMC a C/0 = 4 is typical for non-bipolar PN (IStanghellini et al.ll2009l ). This indicates 
that the models show a higher level of C-enrichment when compared to the observations. 
It is worth mentioning that the model with an enhanced C and O intershell composition, 
computed to reproduce the observational C and O abundances of the M and C-stars in the 
LMC cluster NGC 1978 (refer Section Ej), has a final surface C/0 ratio of C/0 ^ 3.00, in 
good agreement with the LMC PN data. 



Another consequence of the overshoot employed in our models is that the average A 
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Fig. 13. — Mboi versus log T^s for NGC 1978. The black solid-line corresponds to the best-fit 
model. The ellipses mark the predicted C/0 ratios at the end of each TDU episode. A grid of 
the C/0 estimates is shown in the bottom-right corner. The filled black triangles indicate the 
M-stars for which observational C and O abundances exist. The filled green squares indicate 
the C-stars for which observational C and O abundances exist. The size of the ellipse around 
the filled symbols indicates the C/0 ratio of the star. The open black triangles and green 
squares indicate the other M-stars and C-stars, respectively, in the cluster. 
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Fig. 14. — The same as Figure [13] for NGC 1846's best-fit model. 
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values from our full AGB evolution calculations lie in the range 0.40 to 0.70 for our best 
fitting models. The first dredge-up episodes occur for core masses of M™™ ^ 0.56 — O.SSMq. 
Our values can be compared with evolution calculations that do not use any o vershoot at the 



base o f the convective envelope. For direct comparison we use the models of iKarakas et al. 
(I2OO2I ). Figure [15] and Figure [T6] compares our predicted M™™ and Amax (from Table 2) as 
a fun ction of the ini t ial m ass and metallicty to the Z = 0.008 and Z = 0.004 model values 



from iKarakas et al.l (120021 ). We find that, a s expe cted, our models find TDU at smaller 
core masses than the models of IKarakas et al.l (120021 ) without overshoot. Similarly, our A^ax 
values are higher. 

The C-star luminosity function (CSLFs) in the MCs have been used to calibrate the onset 
and efficiency of the TDU in synthetic AGB e volution calculations (iGroenewegen &: de Jong 
I993I : iMarigo et all llOod : llzzard et all 120041 : Istancliffe et aD boosl ). The values derived in 
these four studies are also shown in Figure [15] a nd Figure [T6l The M ™'° values span a 
wide range dow nward from the val ues f ound by IKarakas et al.l (120021 ) and they bracket 
our values. The llzzard et al.l (12004] ) and iMarigo et al.l (Il999l) value seem too smal l to be 
cons istent with our resu lts. In the synthetic calculations by IGroenewegen fc de Jong] (119931 ) 
and Marigo et al. ( 19991), A is a cori s tant p arar neter that is altered to fit the CSLFs but A 
varies in the models of llzzard et al. (2004) and Stancliffe et al. (2005). Given the range of 
masses and metallicities involved in the models and observations, it is difficult to use our 
derived A values to make precise comments on the validity of the A values usined in the 
synthetic models. Observations from clusters with a wider range of initial masses would be 
useful for constraining these parameters as a function of mass and metallicity. 

In Table [2] we give the final core and total masses for the stellar models. Given that we 
expect no further TPs, the final core masses are approximately the final masses, and these can 
be compared to white dwarf masses. Our predicted final masses are ~ 0.59 — 0. 64Mf7>, right at 



the peak of the distribut ion of white dwarf masses in the Milky Way Galaxy (IFerrario et al. 



2OO5I : iLiebert et al.ll2005l ) . Even though these white dwarf masses are derived from a different 
stellar population (notably one more metal rich), there is a strong consistency between our 
predicted masses and the peak of the observed range. 

The main focus of our study is to establish a theoretical explanation for t he observed C, 



O, and F abundances (F abundances only for NGC 1846) in the three clusters. iLederer et al. 



(120091 ) had difficulties finding a theoretical explanation for the NGC 1978 M and C-type 
AGB stars. Our scenario uses a carbon-depleted model with extra-mixing on the RGB along 
with an intershell enhanced in both ^^C and ^^C to explain the composition of the stars 
in NGC 1978. This scenario is somewhat speculative because it is not clear if AGB stars 
experience convective overshoot between the flash-driven convective region and C-0 core. 
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Fig. 15. — The M™™ plotted against initial mass for the Z = .008 (blue dotted- lin e and 
points) and Z = 0.004 (red dashed-line and points) models from ( iKarakas et al.ll2002l ). Also 
plotted are the average M™'° values for the models in each of the three clusters. The values 
of M™™ derived for 1.5Mq models in synthetic AGB calculations designed to reproduce the 
LMC C- star luminosity func tions are shown for lGroenewegen fc de Jond (119931 ) (GD93, open 
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(120041 ) (104, open triangle), and 



(120051 ) (S05, open diamond). 
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Fig. 16. — Same as Fig. [15] but showing the variation of Amax with initial mass for full 
evolution calculations and the values for the synthetic calculations. Note that the A values 
from GD93 and M99 are constant values, while A values vary in 104 and S05 and the value 
plotted is the maximum value. 
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even though this mixing hkely occurs in post-AGB stars to produce the PG 1159 stars. 
For NGC 1846 we find that an a-enhanced initial composition coupled with extra-mixing on 
the RGB followed by standard thermally pulsing AGE evolution reproduces the observed 
abundance patters of the cluster M-stars. We find that extra-mixing , perhaps the beginning 
of hot-bottom burning, is also required on the AGB in order to fit the observational data 
for the C-stars tho u gh we do not model this in our study. This conclusion was also reached 
by iLebzelter et al.l ( l2008l ). For the M-stars in NGC 1846 we insert partial mixing zones 
which turn into the ^^C pocket in our post-processing nucleosynthesis models in order to 
try and re-produce the observed F abundances. We find that a partial mixing zone that is 
a similar size (in mass) to the intershell region at the end of the TP-AGB will reproduce 
the observed F abundance variation with C/0. This result is in contradiction to s-process 
studies that r equire smaller ^^C pocket masses to reproduce the o bserved s-process abundance 
distribution ( IGallino et al.lll998l : iBonacic Marinovic et al.l 120071 ). 

For NGC 419 we have predicted C and O abundances based on a scaled-solar model and 
these results can be compared to the observational abundances once they have been derived. 



DK thanks Dougal Mackey and John Norris for helpful comments on the LMC clusters. 
AIK and DK thank Carlos Iglesias and Forrest Rogers for help with obtaining OPAL tables, 
and Richard Stancliffe for discussions on mixing in giant stars. This work benefited from 
the support of the NCI National Facility at the ANU. 
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Table 1: Parameters used to construct the stellar models for the three target clusters - 
NGC1978, NGC1846, and NGC419. 



Model Type 


MzAMS 


Z 


log Tefr,4 


a 




iVov 




P 


Ti^AGB-tip 




(Mo) 








(Mo) 


(Hp) 




(days) 




NGC1978 


Scaled-Solar 


1.63 


0.006 


3.564 


1.90 


1.56 


2.54 


-4.51 


790 


-5.08 


C-Depleted 


1.63 


0.006 


3.563 


1.90 


1.56 


2.75 


-4.47 


790 


-5.06 


a- Enhanced 


1.63 


0.006 


3.564 


1.90 


1.56 


3.00 


-4.50 


790 


-5.05 










NGCl^ 


m 










Scaled-Solar 


1.86 


0.006 


3.560 


1.74 


1.80 


1.05 


-4.73 


710 


-5.18 


C-Depleted 


1.86 


0.006 


3.558 


1.74 


1.80 


1.05 


-4.77 


710 


-5.18 


a-Enhanced 


1.86 


0.006 


3.559 


1.74 


1.80 


1.41 


-4.78 


710 


-5.17 


NGC419 


Scaled-Solar 


1.90 


0.004 


3.574 


1.74 


1.85 


2.10 


-4.49 


790 


-5.23 



Note: log Tcs^a is the computed log T^e for the M-stars on the AGB at Mboi = —4.00, Nov denotes the 
overshoot parameter in pressure scale heights, Hp (see Section [3.1.2p . M^Jy'^ denotes the computed bolometric 
luminosity at the M/C transition. P is the pulsation period where the superwind mass- loss rate is reached. 
-^bo?^~*'^ denotes the computed AGB-tip bolometric luminosity. 
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Table 2: Characteristic parameters of the AGB models. 



Model Type 


TPs 




Aavg -^dredge 


T^Sell ^c(f) 


Mtot(f) 


Tip(f) 






(Me) 


(Me) 


(lO^K) (Mq) 


(Mo) 


(lO^years) 






NGC1978 


- MzAMS = 1-63, Z = 0.006 






Scaled-Solar 


15 


0.558 0.720 


0.588 0.079 


273 0.605 


0.655 


1.209 


C-Depletcd 


13 


0.555 0.773 


0.653 0.085 


272 0.599 


0.856 


1.398 


a— Enhanced 


13 


0.552 0.817 


0.655 0.097 


272 0.594 


0.647 


1.392 






NGC 1846 


- MzAMS = li 


^6, Z = 0.006 






Scaled-Solar 


16 


0.569 0.551 


0.429 0.051 


279 0.635 


0.933 


1.021 


C-Depleted 


18 


0.573 0.545 


0.433 0.051 


280 0.636 


0.749 


0.978 


a— Enhanced 


17 


0.571 0.609 


0.487 0.060 


279 0.630 


0.705 


1.056 






NGC 419 


- MzAMS = 1-90, Z = 0.004 






Sealed- Solar 


19 


0.555 0.718 


0.608 0.118 


279 0.625 


0.911 


0.997 
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Table 3: Initial, post-FDU, and final C/0 and ^^C/^^C ratios for the three clusters. AGB-tip 
[F/Fe] abundances from the stellar models are also listed. All abundances are by number. 



Model Type 


EM^ 


C/Oi^ 


C/Op-FDu'^ 


C/Ofl^ 12 


C/i3Ci 


12p /13p 


12c/13Cf 


[F/Fe]f 


NGC1978 


Scaled-Solar 


No 


0.501 


0.333 


8.258 


89.4 


23.4 


691.2 


1.147 


C-Depleted 


No 


0.282 


0.185 


7.561 


89.4 


22.7 


1106.2 


1.078 


C-Depleted 


Yes 


0.282 


0.184 


7.601 


89.4 


12.8 


654.0 


1.080 


C-Depleted= 


No 


0.281 


0.184 


7.566 


50.0 


19.3 


948.4 


1.078 


a— Enhanced 


No 


0.282 


0.184 


8.038 


89.4 


22.7 


1175.3 


1.369 


a— Enhanced 


Yes 


0.282 


0.186 


8.079 


89.4 


12.8 


684.4 


1.371 


NGC 1846 


Scaled-Solar 


No 


0.501 


0.309 


4.700 


89.4 


21.9 


386.2 


0.740 


C-Depleted 


No 


0.282 


0.173 


4.480 


89.4 


21.8 


651.5 


0.741 


C-Depleted 


Yes 


0.282 


0.167 


4.489 


89.4 


11.8 


378.6 


0.743 


a— Enhanced 


No 


0.316 


0.196 


4.821 


89.4 


21.8 


621.1 


0.996 


a— Enhanced 


Yes 


0.316 


0.195 


4.832 


89.4 


13.8 


400.8 


0.998 


NGC419 


Scaled-Solar 


No 


0.501 


0.299 


11.954 


89.4 


21.6 


1101.7 


1.393 



"This indicates if extra- mixing is assumed to occur on the RGB. 

''The notations 'i', 'P-FDU', and 'f denote the initial, post-FDU, and final abundance values, respectively. 
"^This represents the carbon-depleted model with an initial ^"^C/^^C = 50 (see Section |4?11 for details). 
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